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Abstract 

Cross sections for e'^p charged current deep inelastic scattering at a centre- 
of-mass energy of 318 GeV have been determined with an integrated lu- 
minosity of 60.9 pb~^ collected with the ZEUS detector at HERA. The 
differential cross sections da/dQ"^, da/dx and da/dy for > 200 GeV^ 
are presented. In addition, d'^a/dxdQ^ has been measured in the kine- 
matic range 280 GeV^ < < 17000 GeV^ and 0.008 < x < 0.42. The pre- 
dictions of the Standard Model agree well with the measured cross sec- 
tions. The mass of the W boson propagator is determined to be 
Mw = 78.9 ± 2.0 (stat.) ± 1.8 (syst.) +^-o (PDF) GeV from a fit to da/dQ^. The 
chiral structure of the Standard Model is also investigated in terms of the (1 — y)^ 
dependence of the double-differential cross section. The structure-function ^2^*-^ 
has been extracted by combining the measurements presented here with previ- 
ous ZEUS results from e~p scattering, extending the measurement obtained in 
a neutrino-nucleus scattering experiment to a significantly higher region. 
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1 Introduction 



Measurements of deep inelastic scattering (DIS) of leptons on nuclei have been vital in 
the development of our understanding of the structure of nucleons. In the Standard 
Model (SM), charged current (CC) DIS is mediated by the exchange of the W boson. 
In contrast to neutral current (NC) interactions, where all quark and antiquark flavours 
participate, only down-type quarks and up-type antiquarks participate at leading order in 
the e~^p CC DIS reaction. Thus, such interactions are a powerful tool for flavour-specific 
investigations of the parton distribution functions (PDFs). Since only left-handed quarks 
and right-handed antiquarks contribute to CC DIS, the distribution of the positron-quark 
centre-of-mass scattering angle, 6*, is a sensitive probe of the chiral structure of the weak 
interaction. 

Measurements of the CC DIS cross sections at HERA have been reported previously by 
the HI [1-4] and ZEUS [5,6] collaborations. These data extended the kinematic region 
covered by fixed-target neutrino-nucleus scattering experiments [7-10] by about two orders 
of magnitude towards larger Q^, the negative square of the four- momentum transfer. 
In addition, the double-differential cross section, d'^a/dxdQ'^, where x is the Bjorken 
scaling variable, was measured for the first time by the HERA collider experiments in e+p 
scattering [11,12] and more recently in e~p scattering [13,14]. 

This paper presents measurements of the e~^p CC DIS single-differential cross sections 
da/dQ"^, da/dx and da/dy, and of d'^a/dxdQ'^ and their comparison with the SM. The 
measurements are based on 60.9 pb~^ of data collected during the running periods in 1999 
and 2000 when HERA collided 27.5 GeV positrons with 920 GeV protons, yielding a centre- 
of-mass energy (-y/s) of 318 GeV. The previous ZEUS measurement [12] at y/s = 300 GeV 
was based on 47.7 pb"'^. The mass of the W boson is determined from the measured cross 
sections and the helicity structure of the Standard Model is investigated. The structure- 
function ^2^^ is extracted and compared to a fixed-target result [15]. 



2 Standard Model prediction 

The CC DIS differential cross section, d^a^^^J dxdQ"^ , for the reaction e"*"p — ^ Vf,X can 
be written at leading order in the electroweak interaction, for longitudinally unpolarised 
beams as [16]: 

^^<.n(e+p) ^ Gl 

dxdQ^ attx (g2 + M^y 

[F+F£%(a;, Q') - y'F^X{x, Q') - Q')]. (1) 
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In this equation, Gp is the Fermi constant, M\y is the mass of the W boson, x is the 
Bjorken scahng variable, y = Q'^/xs (neglecting the masses of the incoming particles) and 
Yj- = 1 ± (1 — y)^. The centre-of-mass energy in the positron-proton collisions is given by 
^/s = 1^ Ef^Ep^ where E^. and E.p are the positron and proton beam energies, respectively. 
The inelasticity, y, is related to 6^* by y = (1 — cos^*)/2. The structure-functions F^'^^^ 
and X-Fg^^p, at leading order in QCD, may be written in terms of sums and differences of 
quark and antiquark PDFs of the proton as follows: 

= AA^^ Q') + Q') + u{x, Q^) + c(x, Q% (2) 

xF^f+p = x[d{x, + g') - g') - c(x, Q% (3) 

where, for example, the PDF d{x, Q"^) gives the number density of down quarks with 
momentum-fraction x at a given . The longitudinal structure function, F^^^^, is zero 
at leading order in QCD. At next-to-leading order (NLO), F^^^^ is non-zero but gives 
a negligible contribution to the cross section except at values of y close to 1, where the 
contribution can be as large as 10%. Since the top-quark mass is large and the off-diagonal 
elements of the CKM matrix are small [17], the contribution from third-generation quarks 
may be ignored [18]. All cross-section calculations presented in this paper were performed 
at NLO in the strong coupling constant. 



3 The ZEUS experiment 

A detailed description of the ZEUS detector can be found elsewhere [19]. A brief outline 
of the components most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [20], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consists of 72 cylindrical drift chamber layers, organised in nine superlayers covering the 
polar-angle^ region 15° < 6 < 164°. The relative transverse- momentum resolution for full- 
length tracks is a{pT)/pT = 0.0058pT© 0.0065 ®0.0014/pt, with pr in GeV. The position 
of the interaction vertex along the beam direction can be reconstructed from the CTD 
tracks with a resolution of about 1 cm in CC events. 

The high-resolution uranium-scintillator calorimeter (CAL) [21] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeter, covering 99.7% 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The polar angle, 9, is measured with respect to the proton beam direction. The 
coordinate origin is at the nominal interaction point. 
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of the solid angle around the nominal interaction point. Each part is subdivided trans- 
versely into towers and longitudinally into one electromagnetic section (EMC) and either 
one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). The smallest 
subdivision of the calorimeter is called a cell. The CAL relative energy resolutions, 
as measured under test-beam conditions, are a{E)/E = 0.18/^/E for electrons and 
a{E)/E = 0.35 /Ve for hadrons, with E in GeV. The timing resolution of the CAL 
is better than 1 ns for energy deposits exceeding 4.5 GeV. The position of the interaction 
vertex along the beam direction can be reconstructed from the arrival time of energy 
deposits in FCAL. The resolution is about 9 cm for events with FCAL energy above 25 
GeV, improving to about 7 cm for FCAL energy above 100 GeV. 

An iron structure that serves as a flux-return path for the magnetic field of the solenoid 
surrounds the CAL and is instrumented as a backing calorimeter (BAG) [22]. Muon 
chambers in the forward, barrel and rear [23] regions are used in this analysis to veto 
background events induced by cosmic-ray or beam-halo muons. 

The luminosity was measured using the Bethe-Heitler reaction ep — > e'yp. The photons 
were detected by the luminosity monitor [24] , a lead-scintillator calorimeter placed in the 
HERA tunnel 107 m from the interaction point in the positron beam direction. Events 
corresponding to an integrated luminosity of 60.9 ± 1.5 pb^^ were used in this analysis. 

4 Monte Carlo simulation 

Monte Carlo (MC) simulations were used to determine the efficiency for selecting events 
and the accuracy of kinematic reconstruction, to estimate the ep background rates and 
to extract cross sections for the full kinematic region. A sufficient number of events was 
generated to ensure that the statistical uncertainties arising from the MC simulation were 
negligible compared to those of the data. The MC samples were normalised to the total 
integrated luminosity of the data. 

Charged current DIS events, including electroweak radiative effects, were simulated us- 
ing the HERACLES 4.6.1 [25] program with the djangoh 1.1 [26] interface to the MC 
generators that provide the hadronisation. Initial-state radiation, vertex and propagator 
corrections and two-boson exchange are included in HERACLES. The mass of the W bo- 
son was calculated using the PDG [17] values for the fine structure constant, the Fermi 
constant, the masses of the Z boson and the top quark, and with the Higgs-boson mass 
set to 100 GeV. The events were generated using the CTEQ5D [27] PDFs. The colour- 
dipole model of ARIADNE 4.10 [28] was used to simulate 0{as) plus leading logarithmic 
corrections to the result of the quark-parton model. As a systematic check, the meps 
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model of lepto 6.5 [29] was used. Both programs use the Lund string model of jetset 
7.4 [30] for the hadronisation. A set of NC events generated with djangoh was used to 
estimate the NC contamination in the CC sample. Photoproduction background was esti- 
mated using events simulated with HERWIG 5.9 [31]. The background from W production 
was estimated using the epvec [32] generator, and the background from production of 
charged-lepton pairs was generated with the lpair [33] and grape [34] programs. 

The ZEUS detector response was simulated with a program based on geant 3.13 [35]. 
The simulated events were subjected to the same trigger requirements as the data, and 
processed by the same reconstruction programs. 

5 Reconstruction of kinematic variables 

The principal signature of CC DIS at HERA is the presence of a large missing transverse 
momentum, PT.miss? arising from the energetic final-state neutrino which escapes detection. 
The quantity PT.miss was calculated from 

^T.miss = Px + Py= Ei sin di cos 0, j + Ei sin Oi sin 0^ 

where the sums run over all CAL energy deposits, Ei (uncorrected in the trigger, but 
corrected [36] for energy loss in inactive material and other effects in the offline analysis), 
and 9i and (pi are the polar and azimuthal angles of the calorimeter deposits as viewed 
from the interaction vertex. The polar angle of the hadronic system, 7/^, is defined by 
COS7. = (P|,^iss - 5')I{Ptmss + 5'), where 5 = Y.E,{1 - cosO,) = j:{E - In the 

i i 

naive quark-parton model, is the angle through which the struck quark is scattered. 
Finally, the total transverse energy, E^, is given by E^ = Yl^i sin6'j. 

i 

The kinematic variables were reconstructed using the Jacquet-Blondel method [37]. The 
estimators of y, and x are: yjB = S/{2Ee), Qjb = -PT,miss/(l " 1/Jb), and xjb = 
<3jb/(^^Jb)- 

The resolution in is about 20%. The resolution in x improves from about 20% at 
X = 0.01 to about 5% at x = 0.5. The resolution in y ranges from about 14% at y = 0.05 
to about 8% at 1/ = 0.83. 

6 Event selection 

Charged current DIS candidates were selected by requiring a large PT,miss and a recon- 
structed event vertex consistent with an ep interaction. The main sources of background 
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come from NC scattering and high-E'j- photoproduction in which the finite energy resolu- 
tion of the CAL or energy that escapes detection can lead to significant measured missing 
transverse momentum. Non-ep events such as beam-gas interactions, beam-halo muons or 
cosmic rays can also cause substantial imbalance in the measured transverse momentum 
and constitute additional sources of background. The selection criteria described below 
were imposed to separate CC events from all backgrounds. 

When the current jet lies in the central region of the detector, i.e. 7/1 is large, tracks 
in the CTD can be used to reconstruct the event vertex, strongly suppressing non-ep 
backgrounds. For CC events with small 7/,, the charged particles of the hadronic final 
state are often outside the acceptance of the CTD. Such events populate the high-x 
region of the kinematic plane. The events were classified first according to 70, the value of 
7/i measured with respect to the nominal interaction point. Subsequently, the kinematic 
quantities were recalculated using the Z-coordinate of the event vertex (Zvtx) determined 
from either CTD tracks, for events with large 70, or the calorimeter-timing information 
for events in which 70 is small. The selection procedures for events with large and small 
7o are described in Sections 6.2 and 6.3, respectively. 

6.1 Trigger selection 

ZEUS has a three-level trigger system [19,38]. At the first level only coarse calorimeter 
and tracking information is available. Events were selected using criteria based on the 
energy, transverse energy and missing transverse momentum measured in the calorimeter. 
Generally, events were triggered with low thresholds on these quantities if a coincidence 
with CTD tracks from the event vertex occurred, while higher thresholds were required 
for events with no CTD tracks. Typical threshold values were 5 GeV (8 GeV) in miss- 
ing transverse momentum, or 11.5 GeV (21 GeV) in transverse energy for events with 
(without) CTD tracks. 

At the second level, timing information from the calorimeter was used to reject events 
inconsistent with the bunch-crossing time. In addition, the topology of the CAL en- 
ergy deposits was used to reject background events. In particular, a tighter cut of 
6 GeV (9 GeV for events without CTD tracks) was made on missing transverse mo- 
mentum, since the resolution in this variable is better at the second level than at the first 
level. 

At the third level, full track reconstruction and vertex finding were performed and used 
to reject candidate events with a vertex inconsistent with an ep interaction. Cuts were 
applied to calorimeter quantities and reconstructed tracks to further reduce beam-gas 
contamination. 
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The efficiency of the full trigger chain for CC DIS events, in the kinematic region of the 
cross-section measurements, was determined using MC simulation to be 96%, and not 
below 78% in any cross-section bin. 

6.2 Offline selection based on a CTD vertex 

Events with 70 > 0.4 rad were required to have a vertex reconstructed from CTD tracks. 
Additional requirements for event selection are described below. 

• Primary vertex: events were required to satisfy |^vtx| < 50 cm. The primary vertex 
reconstructed from the CTD tracks was required to be within the range consistent with 
the ep interaction region. The transverse sizes of the beams were smaller than the CTD 
vertex resolution in X and Y. Furthermore, the average transverse displacement of the 
vertex from its nominal position, as measured by the CTD, was also smaller than the 
resolution. Therefore, the X- and Y-vertex positions were set to zero, their nominal 
positions; 

• Missing transverse momentum: events were required to satisfy PT.miss > 12 GeV and 
-PT.miss > 10 GeV where -Pf.miss is the missing transverse momentum calculated ex- 
cluding the FCAL towers closest to the beam hole. The -Remiss strongly suppresses 
beam-gas events while maintaining high efficiency for CC events; 

• Tracking requirement: tracks associated with the event vertex with transverse mo- 
mentum in excess of 0.2 GeV and a polar angle in the range 15° to 164° were defined 
as "good" tracks. In order to remove beam-gas background, at least one such track 
was required and a cut was also applied in two dimensions on the number of good 
tracks versus the total number of tracks; 

• Rejection of photoproduction: 

— Pt, miss / Et > 0.4 was required for events with Pr.miss < 30 GeV. This cut was 
raised to Pt, miss / Et > 0.55 for events with PT,miss < 20 GeV. This selected events 
with a coUimated energy flow, as expected from a single scattered quark; 

— for events with PT,miss < 20 GeV, the number of good tracks within an azimuthal 
angle of 0.5 rad around the direction of the Pr.miss was counted (A^+) as well as the 
number of tracks opposite to the Pr.miss direction (A^-). The event was rejected if 
the number of tracks in the Pr.miss direction was > 2, or if the asymmetry defined 
as (A^_ - A^+)/(A^_ + A^+) was less than 0.7; 

— for charged current events, there is a correlation between the direction of the 
Pr,miss vector calculated using CTD tracks and that obtained using the GAL. The 
difference between these quantities was required to be less than 0.5 radians for 
Pr.miss < 20 GeV and less than 2.0 radians for Pr.miss > 20 GeV; 
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• Rejection of NC DIS: NC DIS events in which the scattered positron or jet energies 
are poorly measured can have a considerable apparent missing transverse momentum. 
To identify such events, a search for candidate positrons was made using isolated 
electromagnetic clusters in the CAL [39] for events with PT,miss < 30 GeV. Candidate 
positron clusters within the CTD acceptance were required to have an energy above 
4 GeV and a matching track with momentum larger than 25% of the cluster energy. 
Clusters with 9 > 164° were required to have a transverse momentum exceeding 2 GeV. 
Events with a candidate positron satisfying the above criteria and 6 > 30 GeV were 
rejected, since for fully contained NC events, 6 peaks at 2Ee = 55 GeV; 

• Rejection of non-ep background: muon-finding algorithms [40] based on CAL energy 
deposits or muon-chamber signals were used to reject events produced by cosmic rays 
or muons in the beam halo. 

6.3 Offline selection based on a CAL vertex 

For events with 70 < 0.4 rad a vertex reconstructed using calorimeter timing information 
was required. Additional requirements for event selection are described below. 

• Primary vertex: events were required to satisfy |2'vtx| < 50 cm. Zvtx was recon- 
structed from the measured arrival time of energy deposits in FCAL [41]. The X- and 
y-vertex positions were set to zero; 

• Missing transverse momentum: events were required to satisfy Pr.miss > 14 GeV and 
-^Tmiss > 12 GeV. The conditions on missing transverse momentum were tightened, 
compensating for the relaxation of the track requirements; 

• Rejection of non-ep background: the muon-rejection cuts described in Section 6.2 
were used. A class of background events arose from beam-halo muons that produced a 
shower inside the FCAL. To reduce this background, topological cuts on the transverse 
and longitudinal shower shape were imposed; these cuts rejected events in which the 
energy deposits were more coUimated than for typical hadronic jets. 

6.4 Final event sample 

To restrict the sample to kinematic regions where the resolution in the kinematic variables 
was good and the backgrounds small, the requirements Qjg > 200 GeV^ and ?/jb < 0.9 
were imposed. 

All events were visually inspected, and five cosmic-ray and halo-muon events were removed 
from the large-70 sample and 11 from the small-70 sample. A total of 1164 data events 
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satisfied these criteria in the large-70 sample and 292 in the small-70 sample, in good 
agreement with the 1183 and 285 predicted by the MC simulation for the large- and 
small-70 samples, respectively. 

Figure 1 compares the distributions of data events entering the final CC sample with 
the MC expectation for the sum of the CC signal and ep background events. The MC 
simulations give a good description of the data. 

7 Cross section determination and systematic 
uncertainties 



The measured cross section in a particular kinematic bin, for example for d'^a/dxdQ'^, was 
determined from 

d^a TVdata-iVbg 



dxdQ"^ Nmc dxdQ"^ 

where A^'data is the number of data events, A^'bg is the number of background events es- 
timated from the MC simulation and A^mc is the number of signal MC events. ji^ncdQ'^ 
is the Standard Model prediction evaluated in the on-shell scheme [17] using the PDG 
values for the electroweak parameters and the CTEQ5D PDFs [27]. A similar procedure 
was used for da/dQ"^, da/dx and da/dy. Consequently, the acceptance, as well as the 
bin-centring and radiative corrections were all taken from the MC simulation. The cross 
sections da /dQ'^ and da/dx were extrapolated to the full y range using the SM predictions 
calculated with the CTEQ5D PDFs. The extrapolation factors ranged from 3% at the 
lowest to 28% at the highest in da jdQ^ and from zero to 8% in da/dx. 

The systematic uncertainties in the measured cross sections were determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of the 
cross sections. 

• Uncertainties of the calorimeter energy scale: 

— the relative uncertainty of the hadronic energy scale was 2% for the RCAL and 1% 
for the FCAL and BCAL [14]. Varying the energy scale of the calorimeter sections 
by these amounts in the detector simulation induces small shifts of the Jacquet- 
Blondel estimators of the kinematic variables. The variation of the energy scale 
of each of the calorimeters simultaneously up or down by these amounts gave the 
systematic uncertainty on the total measured energy in the calorimeter. This was 
found to give shifts in the cross sections which were correlated between kinematic 
bins in Tables 1 and 2); 
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— the uncertainty in the cross sections from the effect of the energy scale on the 
measurement of 7?,, was obtained by increasing (decreasing) the FCAL and RCAL 
energy scales together while the BCAL energy scale was decreased (increased). 
This gave shifts in the cross sections that were correlated between kinematic bins 
(^2 in Tables 1 and 2); 

— the sensitivity of the cross-section measurements to the fractions of the energy 
deposited in the EMC and HAC sections of the calorimeter was determined by si- 
multaneously increasing the energy measured in the EMC section of the calorime- 
ter by 2% and decreasing the energy measured in the HAC section by 2%, and 
vice-versa. This was done separately for each of the calorimeter sections; 

The final systematic error attributed to the uncertainty in the hadronic energy scale 
was obtained by taking the quadratic sum of these three types of estimate. The 
resulting systematic shifts in the measured cross sections were typically within ±5%, 
but increased to ±20% in the highest and x bins; 

• Energy leakage: of the accepted events, 4% have a measurable energy leakage from 
the CAL into the BAC. The average leakage of transverse energy for these events is 
5% of that observed in the CAL. Both the fraction of events with leakage and the 
average amount of leakage are well modelled by the MC simulation and the effect on 
the cross-section measurement is negligible; 

• Variation of selection thresholds: the selection thresholds were varied by the typical 
resolution of the CAL quantities (±10%) in both data and MC and no significant 
systematic effects were observed. The criteria for good tracks was also varied and no 
systematic effects in the measured cross sections were observed; 

• Parton-shower scheme: the MEPS model of LEPTO was used instead of the ARIADNE 
model. This gave shifts in the cross sections which were found to be correlated between 
kinematic bins (^3 in Tables 1 and 2). The largest effects were observed in the bins at 
high (±20%) and highest and lowest x (±6%). The largest effect in the double- 
differential cross section was seen in the low-Q^, low-x bins, where it amounted to 
±10%; 

• Background subtraction: the uncertainty in the small contribution from photoproduc- 
tion was estimated by fitting a linear combination of the Pt^ss/ distributions of the 
signal and the background MC samples to the corresponding distribution in the data, 
allowing the normalisation (iVphp) of the photoproduction MC events to vary. No cut 
on PT^miss/ Et was applied for this check. Varying A^php by ±20%, corresponding to 
twice the uncertainty given by the fit, results in modifications of the cross sections 
within ±2%; 

• Parton distributions: the CC MC events were generated using the CTEQ5D PDFs [27]. 
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The ZEUS-S fit [42] was used to examine the influence of variations of the PDFs on 
the cross-section measurement through differences in the acceptance and bin-centring 
corrections. The Monte Carlo events were reweighted to the extremes of the cross- 
section prediction allowed by the fit. The change in the measured cross section was 
typically < 1%, except at high where it was —5% and at high x where it was +4%; 

• NLO QCD corrections: the DJANGOH program neglects the Fl contribution and NLO 
QCD corrections to X-F3 when generating CC events. The corresponding effect on the 
cross-section measurement was evaluated by reweighting the MC events with the ratio 
of the cross sections with and without NLO QCD corrections. The largest effect, of 
—3%, was observed in the highest y bin. 

The uncertainties associated with the trigger and the measurement of the vertex positions 
were negligible. The individual uncertainties were added in quadrature separately for the 
positive and negative deviations from the nominal cross-section values to obtain the total 
systematic uncertainties listed in Tables 1 and 2. The 0{oi) electroweak corrections to CC 
DIS have been discussed by several authors [43,44]. Various theoretical approximations 
and computer codes gave differences in the CC cross sections of typically ±(1 — 2)% 
or less. However, the differences can be as large as ±(3 — 8)% at high x and high y. 
This uncertainty and the uncertainty of 2.5% on the measured total luminosity were not 
included in the figures and the tables of the cross sections. 

8 Results 

The total cross section for e^'p CC DIS in the kinematic region > 200 GeV^ is 

(^to?(Q^ > 200 GeV^) = 34.8 ± 0.9(stat.)i?;^(syst.) pb. 

In this case the uncertainty in the measured luminosity is included in the systematic 
uncertainty. The result is in agreement with the SM expectation of 37.01;q g pb evaluated 
using the ZEUS-S fit. 

8.1 Single-differential cross sections 

The single-differential cross sections da/dQ"^^ da/dx and da/dy for > 200 GeV^ are 
shown in Figs. 2, 3 and 4, respectively, and compiled in Table 1 including details of 
the systematic uncertainties that are correlated between cross-section bins. The SM 
cross sections derived from Eq. (1) using the ZEUS-S fit, the CTEQ6D [45] and the 
MRST(2001) [46] parameterisations of the PDFs are shown, together with the ratios of 
the measured cross sections to the SM cross section evaluated with the ZEUS-S fit. 
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The cross sections do jdQ^ and da /dx drop by four and three orders of magnitude, respec- 
tively, due to the effect of the W^-boson propagator and the decreasing quark density at 
large x. The ZEUS-S fit was based on fixed-target DIS data obtained at much lower 
(< 100 GeV^) and from ZEUS NC data at large Q^. The excellent description of the data 
by the SM prediction based on this fit confirms both the decomposition of the proton mo- 
mentum into different quark flavours, specifically the down-quark contribution, and the 
evolution of parton distributions towards scales considerably larger than the py-boson 
mass. At very large x and Q^, the uncertainty in the prediction derived from the ZEUS-S 
fit, and also the global fits, reflects the lack of data constraining the ci-quark density. 

8.2 Double-differential cross sections 

The reduced double-differential cross section, a, is defined by 



At leading order in QCD, a{e^p — ^ VeX) depends on the quark momentum distributions 
as follows: 



The reduced cross sections are displayed as functions of and x in Figs. 5 and 6, 
respectively, and compiled in Table 2 including details of the systematic uncertainties 
that are correlated between cross-section bins. The predictions of Eq. (1), evaluated 
using the ZEUS-S fit, the CTEQ6D and the MRST(2001) PDFs give a good description 
of the data. The contributions from the PDF combinations {d + s) and {u + c), obtained 
in the MS scheme from the ZEUS-S fit, are shown separately in Fig. 6. 

8.3 Helicity studies 

The W boson couples only to left-handed fermions and right-handed antifermions. There- 
fore, the angular distribution of the scattered quark in e^q CC DIS will be fiat in 6'*, while 
it will exhibit a (1 + cos 9*Y distribution in e^q scattering. Since (1 — y)^ oc (1 + cos 6'*)^, 
the helicity structure of CC interactions can be illustrated by plotting the reduced double- 
differential cross section of Eq. (4) versus {1—yY in bins of x. In the region of approximate 
scaling, i.e. x ~ 0.1, this yields a straight line. At leading order in QCD, the intercept of 
this line gives the {u + c) contribution, while the slope gives the {d + s) contribution. 

Figure 7 shows a as a function of (1 — yY for the e^p CC DIS data, compared to the 
previously published e~p data [14]. At large x, e~^p CC DIS is sensitive to the valence part 




a{e^p VeX) = x[u + c+{l- yf{d + s)] . 



(4) 
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of d{x,Q'^), while e p CC DIS is sensitive to the valence part of u{x,Q'^). Equivalently 
to Eq. (4), the reduced cross section for e~p CC DIS can be written as 



permitting a similar interpretation for the intercept and slope in terms of the appropriate 
parton densities. Scaling violations can be observed in the theoretical prediction as (1— y)^ 
approaches one. The data agree with the expectation of the SM from the ZEUS-S fit. 

8.4 Mass of the W boson 

The fall in the cross section da/dQ"^ with increasing depends on M^/{Cp + M^)^. 
Fitting da/dQ"^ with Gp fixed at the PDG [17] value of 1.16639 ■ 10"^ GeV'^ using the 
ZEUS-S fit PDFs and My/ treated as a free parameter, gives: 



where the third uncertainty was estimated by varying the PDFs within the uncertain- 
ties given by the ZEUS-S fit. The systematic uncertainty includes contributions from 
the sources identified in Section 7 and the uncertainty on the measured luminosity. This 
measurement, in the space-like region, is in good agreement with the more precise mea- 
surements of W-hoson mass in the time-like region [17]. 



It is possible to extract the quark singlet distribution, ^2^*", by combining the measure- 
ments presented in this paper with the previous ZEUS e~p CC DIS results [14]. This can 
be compared to the corresponding result from fixed-target neutrino-nuleon scattering. 
Equations 2 and 3 show that only down-type quarks and up-type antiquarks contribute 
to e^p CC DIS in the Standard Model. The cross section for e~p CC DIS is given by 



cr(e p VeX) = X [m -|- c -|- (1 — yY{d + s)] 



Mw = 78.9 ± 2.0(stat.) ± 1.8(syst.):f g (PDF) GeV, 



8.5 Extraction of 





x[u{x, + c(x, + d{x, + s(x, Q% 



xFl 



.CC 
3,e~p 
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Therefore, the sum of and F^^^p represents the contribution from all quark and 

antiquark flavours and can be extracted from the measured CC cross sections as 

' ~ Gl\ Ml ) y\ dxdQ^ + dxdQ^ J+iX[x^3,^L), 

where A(xF3, Fl) denotes a correction term taking into account the xF^ and Fl structure 
functions. The correction is given by 

The dominant uncertainty on the extracted structure-function is statistical. The system- 
atic uncertainties 6i, 62, S3 (defined in Section 7) and the luminosity uncertainties were 
considered fully correlated between the two data sets. The other systematic uncertainties 
were treated as uncorrelated. Figure 8 shows the extracted structure-function Fj^*-^ as a 
function of for different values of x. The size of the correction A^xF^, Fl), computed 
at NLO in QCD using the ZEUS-S fit, is shown as a shaded area in each bin. It can 
be seen that the correction term is smaller than the uncertainties on the measurement 
for low values of x but becomes sizeable at higher values of x and Q^. The uncertainity 
on the correction A(xF3,Fl) was also computed using the ZEUS-S fit and found to be 
small. The corresponding result for z/Fe interactions from the CCFR Collaboration [15] 
is shown, after correcting for heavy-target effects [47]. It can be seen that both sets of 
results, spanning more than four orders of magnitude in Q^, are well described by the 
SM prediction evaluated using the ZEUS-S fit. The prediction is somewhat above the 
CCFR data in the highest x bins. Note that the CCFR F^'^ results were not included in 
the ZEUS-S fit. The global fits MRST(2001) and CTEQ6D which use the CCFR results 
obtain a reasonable description of both data sets over the entire kinematic range. 

9 Summary 

Differential cross sections for charged current deep inelastic scattering, e'^p — > Pe^, have 
been measured for > 200 GeV^ using 60.9 pb~^ of data collected with the ZEUS 
detector during the period 1999 to 2000. The double-differential cross section d'^a/dx dQ^ 
is presented in the kinematic range 280 GeV^ < < 17 000 GeV^ and 0.008 <x < 0.42. 
The chiral structure of the Standard Model was investigated by plotting the double- 
differential cross section as a function of (1 — y)^. The mass of the W boson, determined 
from a fit to da/dQ^, is Mw = 78.9 ± 2.0 (stat.) ± 1.8 (syst.) (PDF) GeV. The singlet 
structure-function Fg^*" has been extracted for the first time at HERA, by combining the 
measurements presented here with previous ZEUS measurements. The Standard Model 
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gives an excellent description of all data, confirming the decomposition of the proton 
momentum into different quark flavours, specifically the down-quark distribution, and 
also verifies the QCD evolution of parton distributions towards scales considerably larger 
than the H^-boson mass. 
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range 
(GcV2) 



(GeV2) 



JVda 



bg 



da/dQ'^ 

<5unc(%) I 



52(%) 



<53{%) 



da/dQ'^ (pb/GeV^) 



0.13x 
0.15-' 



200 ■ 
400 ■ 
711 - 
1265 ■ 
2249 ■ 
4000 ■ 
7113 - 
12649 ■ 
22494 ■ 



400 
711 
1265 
2249 
4000 
7113 
12649 
22494 
60000 



280 

530 

950 

1700 

3000 

5300 

9500 

17000 

30000 



159 
204 
306 
324 
235 
155 
59 
11 
3 



14.2 
2.3 
4.2 
1.4 
0.9 
0.2 
0.2 
0.0 
0.0 



+ 2.0 
-3.1 
+ 1.8 
-1.5 
+ 1.0 
-1.2 
+ 1.2 
-0.9 
+ 1.1 
-1.0 
+0.8 
-1.1 
+ 1.8 
-1.1 
+ 1.5 
-2.7 
+4.7 
-6.3 



+ 1.9 

-1.8 

+ 1.7 

-0.8 

+0.7 

-1.0 

+0.5 

-0.0 

-1.9 

+ 1.8 

-2.2 

+2.9 

-6.4 

+5.9 

-8.8 

+ 11 

-16 

+ 18 



+0.1 
-0.5 
+0.8 
+0.8 
-0.5 
-0.5 
+0.3 
+0.2 
-0.1 
+0.2 
+0.1 
+0.6 
+0.3 
-1.9 
+2.5 
-2.4 
+5.1 
-7.8 



-3.8 
+3.8 
-1.6 
+ 1.6 
-2.7 
+2.7 
+ 1.6 
-1.6 
-0.6 
+0.6 
+0.4 
-0.4 
+2.6 
-2.6 
+6.5 
-6.5 
+21 
-21 



(2.85 ±0.22: 
(1.81 ±0.13^0 
(1.30 ±0.08^0 
(7.16±0.4llg 



(2.90 ±0.19] 
(1.07 ±0.09^0 
(2.20±0.29lo 
O nc;+0-82 +0 

fr, -19 + 2.06 +0 

„o 



06 
1.04/ 
1.04-1 
.04/ 
.15-, 
.13/ 
1.06-1 
).07/ 
).03-i 
1.03/ 
1.15-1 
).16/ 

;f3).10-5 



10-2 
10-2 
10-2 
10-3 
10-3 
10-3 
10-"* 



da/dx 



X ranee 



JVda 



bg 



c(%) 



<52(%) 



<53{%) 



da/dx (pb) 



0.010 -0.022 
0.022 -0.046 
0.046 -0.100 
0.100 -0.178 
0.178 -0.316 
0.316 -0.562 
0.562 -1.000 



0.015 
0.032 
0.068 
0.130 
0.240 
0.420 
0.650 



167 
351 
425 
258 
173 
45 
2 



7.7 
6.2 
4.8 
0.9 
0.5 
0.1 
0.0 



+ 1.0 
-1.3 
+ 0.5 
-0.7 
+0.6 
-0.9 
+0.6 
-0.6 
+0.8 
-0.5 
+ 1.1 
-0.8 
+6.5 
-2.1 



+2.2 
-0.9 
+0.4 
-0.4 
-0.5 
+0.6 
-0.8 
+ 1.0 
-2.2 
+2.4 
-4.3 
+4.1 
-12 
+ 14 



+0.1 
-0.0 
+0.1 
+0.2 
+0.7 
-0.6 
+0.9 
-0.6 
-1.3 
+ 1.3 
-3.2 
+3.1 
-11 
+ 16 



-5.5 
+5.5 
+0.9 
-0.9 
+0.3 
-0.3 
+0.9 
-0.9 
+0.7 
-0.7 
-3.7 
+3.7 
-8.7 
+8.7 



(4.58±0.36t»;2«) 

{2.92±0.Wto°l) 

(1.59 ±0.081^2) 

(7.22±OA5to\l) 
+o!o9-i 

0.08^ 

(5.98 ±0.90lo 39) 

C4 43+5.84 +1.06-, . 
('i-'iO_2.S7 -0.81^ 



(3.01 ±0.23 



■ 102 

■ 102 

■ 102 

•lOi 
■lOi 

■ 10" 

10-1 



dcr/dy 



y range 



y 



bg 



c(%) 



Sl(%) 



S2{%) 



<53{%) 



dcr/dy (pb) 



0.00 ■ 
0.10 ■ 
0.20 ■ 
0.34 ■ 
0.48 ■ 
0.62 ■ 
0.76 ■ 



0.10 
0.20 
0.34 
0.48 
0.62 
0.76 
0.90 



0.05 
0.15 
0.27 
0.41 
0.55 
0.69 
0.83 



264 
360 
316 
219 
146 
102 
49 



7.0 
6.3 
4.1 
1.9 
2.3 
1.4 
0.4 



+0.8 
-1.2 
+0.5 
-0.8 
+0.5 
-0.8 
+0.7 
-0.6 
+ 1.0 
-1.0 
+ 1.5 
-1.5 
+3.0 
-3.0 



+0.5 
-0.6 
+0.4 
-0.2 
-0.3 
+0.2 
-0.0 
+0.5 
-0.2 
+0.7 
-2.3 
+3.1 
-4.2 
+4.3 



+0.2 
-0.1 
-0.3 
+0.4 
-0.8 
+0.7 
-0.1 
+0.5 
+ 1.0 
-0.5 
+2.6 
-2.5 
+2.9 
-2.4 



+ 1.7 
-1.7 
-1.0 
+ 1.0 
-3.5 
+3.5 
-0.2 
+0.2 
-2.2 
+2.2 
-1.8 
+ 1.8 
-1.8 
+ 1-8 



(7.79±0.48lo;i7) 



(6.86 ±0.37 



+0.08-1 
0.09^ 

(4.46±0.25l0;i«) 

(3.36±0.23lgg3) 

+0.07 -, 
0.06^ 



(2.58 ±0.21 
(2.15±0.22lO;J^) 



(1.53 ±0.22: 



10^ 
lOi 
lOi 
lOi 
lOi 
10^ 
lOi 



Table 1: Values of the differential cross sections da/dQ"^, dcr/dx and da/dy. 
The following quantities are given: the kinematic range of the measurement; the 
value of x,y and at which the cross section is quoted; the number of data events; 
the number of expected background events; the uncorrelated systematic uncertainty; 
those systematic uncertainties with correlations between cross-section bins 5i, 82 
and (^3 defined in Section 7 and the measured cross section, with statistical and total 
systematic uncertainties. The uncertainty on the measured luminosity of 2.5% is 
not included in the total systematic uncertainty. 
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Q'- ( GeV ) 


X 




^bg 


l5unc(%) 


<5i (%) 


<52(%) 




(J 


280 





008 


26 


2.5 


+2.8 
-2.8 


+4.4 
-2.5 


+0.3 
+2.1 


— 5.3 
+5.3 


(1.47±0.28+'^;JJ) 


■ 10° 


280 





015 


49 


5.3 


+ 1.9 
-3.0 


+ 1.4 
-2.2 


— 0.6 
-0.1 


— 10 
+ 10 


(1.08±0.15l°;i^) 


■ 10° 


280 





032 


55 


3.3 


+ 1.5 
-1.9 


+2.1 
+0.2 


+ 1.6 

-0.6 


—3 1 
+3.1 


(8.52±1.12tO-|I) 


10" 


— 1 


280 





068 


24 


2.8 


+0.8 
-5.1 


+0.6 
-4.0 


-0.7 
-1.7 


— 1.0 

+ 1.0 


(4.14±0.78t°;^^) 


10" 




530 





015 


52 


0.5 


+ 1.0 
-1.0 


+2.1 
-0.7 


— 0.9 
+0.3 


— 7.6 

+7.6 


(8.39±1.18t°j^) 


10" 




530 





032 


57 


0.7 


+ 1.0 
-0^8 


+0.5 

-1.3 


— 0.2 
+0.9 


+3.7 
-3^7 


(6.15±0.82tO-J4) 


10" 




530 





068 


59 


0.6 


+ 1.8 
-1.3 


+ 1.6 
+0.5 


+2.1 
+0.9 


— 2.4 

+2.4 


(6.28 ± 0.83^0.??) 


10" 




530 





130 


25 


0.2 


+7.2 
-1.1 


+4.2 
+ 1.2 


+4.1 
+2.5 


+10 
-10 


(4.40±0.89tg;°^) 


10" 




950 





015 


52 


2.0 


+2.0 
-1.6 


+ 1.1 
+0.6 


+ 1.5 

-0.9 


— 8.5 

+8.5 


(6.85 ± 0.94t° Jg) 


10" 




950 





032 


102 


0.8 


+0.5 
-0.6 


+ 1.5 

-0.7 


-0.4 
+0.5 


— 5.4 

+5.4 


(6.97±0.70to3g) 


10" 




950 





068 


84 


0.8 


+0.5 
-1.6 


— 0.4 
-0.9 


— 1.0 
-0.2 


+2.0 
-2.0 


(5.63±0.62Io:tg) 


10" 


_ 1 


950 





130 


48 


0.2 


+0.5 
-4.0 


+0.6 
-3.8 


— 0.4 
-2.7 


+ 1.4 
-1.4 


(4.86±0.71Io3i) 


10" 


-1 


950 





240 


20 


0.3 


+5.1 
-1.1 


+0.9 
-0.7 


— 4.5 
-0.4 


—4.4 
+4.4 


(2.87 ±0.64^0.19) 


10" 


-1 


1700 





032 


105 


0.9 


+0.9 
-0.8 


-1.0 
-0.6 


+0.4 
-0.0 


+3.2 
-3.2 


(5.55±0.55to.l9) 


10" 


-1 


1700 





068 


105 


0.2 


+ 1.2 
-0.5 


+ 1.1 
+0.4 


-0.1 
+0.5 


-1.0 
+ 1.0 


(4.84±0.48toJg) 


10" 


-1 


1700 





130 


57 


0.2 


+ 1.3 
-0.6 


+0.4 
+0.9 


+ 1.5 

-0.2 


+3.3 
-3.3 


(3.71±0.50t°;i3) 


10" 


-1 


1700 





240 


39 


0.1 


+ 1.3 
-3.4 


-0.4 
+0.4 


-0.5 
+0.2 


-3.6 
+3.6 


(2.73±0.44+°-;t°) 


10" 


-1 


3000 





068 


97 


0.3 


+0.9 
-1.6 


-2.2 
+ 2.2 


+2.1 
-1.3 


-1.1 

+ 1.1 


(3.71±0.39tg;t^) 


10" 


-1 


3000 





130 


55 


0.1 


+ 1.0 
-0.5 


-2.0 
+ 1.6 


-1.4 
+ 1.9 


-3.2 
+3.2 


(2.73±0.37toii) 


10" 


-1 


3000 





240 


44 


0.1 


+2.3 
-0.5 


-2.1 
+2.1 


-1.3 
+1.8 


-3.0 

+3.0 


(2.17±0.33tO;J«) 


10" 


-1 


3000 





420 


7 


0.0 


+3.5 
-1.3 


-0.4 
+ 1.4 


-2.3 
+ 1.7 


-3.1 
+3.1 


r,.+2.29 +0.22 N 
V^-^^-1.57 -0.17/ 


10- 


2 


5300 





068 


49 


0.1 


+ 1.0 
-1.4 


-2.7 
+2.3 


+1.3 
-2.5 


+2.3 
-2.3 


(2.28±0.33t°;?g) 


10" 


-1 


5300 





130 


45 


0.0 


+0.8 
-0.8 


-1.4 
+2.5 


+ 1.0 
+0.1 


-1.2 
+ 1.2 


(1.96±0.29t^;!54) 


10" 


-1 


5300 





240 


41 


0.0 


+0.9 
-0.7 


-2.4 
+3.7 


-2.0 
+4.4 


-0.7 
+0.7 


(1.80±0.28t°;Ji) 


10" 


-1 


5300 





420 


20 


0.0 


+ 1.1 
-3.9 


-2.7 
+3.8 


-2.4 
+4.2 


+3.6 
-3.6 


(9.85±2.22t°-g) 


10" 


-2 


9500 





130 


21 


0.1 


+3.5 
-2.4 


-6.7 
+5.8 


+5.2 
-5.7 


+3.9 
-3.9 


(9.99±2.18tO|*) 


10" 


-2 


9500 





240 


22 


0.0 


+ 1.0 
-2.2 


-4.2 
+3.6 


-3.5 
+1.3 


+3.6 
-3.6 


(1.03 ±0.22+° °^) 


10" 


-1 


9500 





420 


8 


0.0 


+3.4 
-0.6 


-7.6 
+6.1 


-6.2 
+5.7 


+0.3 
-0.3 


1% 7S+1.86 +0.34 
'°-1.31 -0.37.* ■ 


10" 


2 


17000 





240 


3 


0.0 


+3.9 
-2.8 


-8.8 
+ 11 


+6.1 
-5.0 


+6.5 
-6.5 


n c:c + 1.49 +0.23-, 
ll-OO_0.84 -0.19^ ■ 


10" 


2 


17000 





420 


6 


0.0 


+0.7 
-3.8 


-7.0 
+7.8 


-6.5 
+5.6 


+2.8 
-2.8 


/o in+l-91 +0.32-, 
VO.lM^j^ 27 -0.34^ ■ 


10- 


2 



Table 2: Values of the reduced cross section. The following quantities are given: 
the values of and x at which the cross section is quoted; the number of data 
events; the number of expected background events; the uncorrelated systematic un- 
certainty; those systematic uncertainties with correlations between cross-section 
bins 5i, 82 and 5^ defined in Section 7 and the measured cross section, with statisti- 
cal and total systematic uncertainties. The uncertainty on the measured luminosity 
0/2.5% is not included in the total systematic uncertainty. 
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Figure 1: Comparison of the final e^p CC data sample (solid points) with the 
sums of the signal and ep-background Monte Carlo simulations (light shaded his- 
tograms). The simulated ep-background events are shown as the dark shaded his- 
tograms. Shown are the distributions of (a) the missing transverse momentum, 
-Pr.miss; (b) -Pr.miss excluding the very forward cells, P-rmiss' (^) variable S, de- 
fined in Section 5, (d) the ratio of missing transverse momentum to total transverse 
energy, PT,miss/ Et, (g) ■jh, (f) the number of good tracks, (g) the Z position of the 
CTD vertex for the large-'-fQ sample and (h) the Z position of the timing vertex for 
the small-'jo sample. 
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Figure 2: (a) The e^p CC DIS Born cross section da/dQ"^ for data and the Stan- 
dard Model expectation evaluated using the ZEUS-S, the CTEQ6D and the MRST 
(2001) PDFs. The data are shown as the filled points, the statistical uncertanties 
are indicated by the inner error bars (delimited by horizontal lines) and the full er- 
ror bars show the total uncertainty obtained by adding the statistical and systematic 
contributions in quadrature, (b) The ratio of the measured cross section, da/dQ"^, 
to the Standard Model expectation evaluated using the ZEUS-S fit. The shaded band 
shows the uncertainties from the ZEUS-S fit. 
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Figure 3: (a) The e~^p CC DIS Born cross section da/dx for data and the Standard 
Model expectation evaluated using the ZEUS-S, the CTEQ6D and the MRST (2001) 
PDFs. The data are shown as the filled points, the statistical uncertanties are 
indicated by the inner error bars (delimited by horizontal lines) and the full error 
bars show the total uncertainty obtained by adding the statistical and systematic 
contributions in quadrature, (b) The ratio of the measured cross section, da/dx, to 
the Standard Model expectation evaluated using the ZEUS-S fit. The shaded band 
shows the uncertainties from the ZEUS-S fit. 
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Figure 4: (a) The e^p CC DIS Born cross section dcr/dy for data and the Standard 
Model expectation evaluated using the ZEUS-S, the CTEQ6D and the MRST (2001) 
PDFs. The data are shown as the filled points, the statistical uncertanties are 
indicated by the inner error bars (delimited by horizontal lines) and the full error 
bars show the total uncertainty obtained by adding the statistical and systematic 
contributions in quadrature, (b) The ratio of the measured cross section, da/dy, to 
the Standard Model expectation evaluated using the ZEUS-S fit. The shaded band 
shows the uncertainties from the ZEUS-S fit. 
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Figure 5: The reduced cross section, a, as a function of , for different fixed 
values of X. The data are shown as the filled points, the statistical uncertanties are 
indicated by the inner error bars (delimited by horizontal lines) and the full error 
bars show the total uncertainty obtained by adding the statistical and systematic 
contributions in quadrature. The expectation of the Standard Model evaluated using 
the ZEUS-S, the CTEQ6D and the MRST(2001) PDFs is shown by the solid, dashed 
and dotted lines, respectively. The shaded band shows the uncertainties from the 
ZEUS-S fit. 
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Figure 6: The reduced cross section, a, as a function ofx, for different values of 
. The data are shown as the filled points, the statistical uncertanties are indicated 
by the inner error bars ( delimited by horizontal lines ) and the full error bars show 
the total uncertainty obtained by adding the statistical and systematic contributions 
in quadrature. The expectation of the Standard Model evaluated using the ZEUS- 
S fit is shown as a solid line. The shaded band shows the uncertainties from the 
ZEUS-S fit. The separate contributions of the PDF combinations x{u + c) and 
(1 — yYx{d + s) are shown by the dotted and dashed lines, respectively. 
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Figure 7: The reduced cross section, a, as a function of (1 — y)^, for different 
fixed values of x, for e'^p (solid points) and e~p (open circles) CC DIS. The data 
are shown as the points, the statistical uncertanties are indicated by the inner error 
bars (delimited by horizontal lines) and the full error bars show the total uncertainty 
obtained by adding the statistical and systematic contributions in quadrature. The 
expectation of the Standard Model evaluated using the ZEUS-S fit is shown as a 
solid line. The contributions of the PDF combinations x{u + c) and x{u + c) are 
shown by the dashed and dotted lines, respectively. 
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Figure 8: The structure-function as a function ofQ^, for different fixed values 
ofx, extracted from ZEUS data (solid points) and compared to CCFR measurements 
corrected for heavy-target effects (open circles). The expectation of the Standard 
Model evaluated using the ZEUS-S fit is shown as a solid line. The shaded hand 
shows the size of the correction AlxF^, Fi,) defined in Section 8.5. 



27 



